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Spectral breaks as a signature of cosmic ray induced turbulence in the Galaxy 
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We show that the complex shape of the cosmic ray (CR) spectrum, as recently measured by 
PAMELA and inferred from Fermi-LAT 7-ray observations of molecular clouds in the Gould belt, 
can be naturally understood in terms of basic plasma astrophysics phenomena. A break from a 
harder to a softer spectrum at blue rigidity _R ~ 10 GV follows from a transition from transport 
dominated by advection of particles with Alfven waves to a regime where diffusion in the turbulence 
generated by the same CRs is dominant. A second break at i? ~ 200 GV happens when the diffusive 
propagation is no longer determined by the self-generated turbulence, but rather by the cascading 
of externally generated turbulence (for instance due to supernova (SN) bubbles) from large spatial 
scales to smaller scales where CRs can resonate. Implications of this scenario for the cosmic ray 
spectrum, grammage and anisotropy are discussed. 
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Introduction — The spectrum of cosmic rays (CRs) as 
observed at the Earth is not a perfect power law, even 
below the knee, at ~ 10^^ eV. Recent measurements by 
the PAMELA experiment have shown that the spectrum 
has a change of slope at blue ~ 230 GV, e.g., for pro- 
tons, from oc £'-2-85 for E < 230 GeV to blue oc E'^-^'^ 
for E > 230 GeV [J. Moreover, Fermi-LAT 7-ray data 
between 100 MeV and 100 GcV from molecular clouds in 
the Gould belt (but off the Galactic disc) show evidence 
for a composite CR spectrum, with a slope a « 1.9 be- 
low - 10 GeV and a « 2.9 between 10 and 200 GeV 
Q. This analysis has recently been repeated by 3 who 
find very similar results. The spectrum of CR protons 
as measured by PAMELA in the highest energy bins is 
compatible with balloon measurements by CREAM 
At energies below a few tens of GeV the spectrum ob- 
served by PAMELA is affected by solar modulation and 
hence it is non-trivial to use these data for comparison 
with the findings of 

These spectral features are so far either attributed to 
propagation effects from nearby sources or modeled with 
corresponding breaks in either the injection spectrum of 
CRs, or in the diffusion coefficient experienced by these 
particles on their way through the Galaxy Since the 
CR fiux at Earth is dominated by distant sources, it ap- 
pears rather unlikely that nearby sources may produce 
fluctuations of order unity at these energies 0]. Recently 
proposed that the change of spectrum may be induced 
by a diffusion coefficient that is not separable in energy 
and space. 

In this Letter we propose that the complex observed 
CR spectrum is telling us something about how and why 
particles diffuse in the interstellar medium (ISM) and de- 
rives from diffusion of CRs on a background of waves 
partly due to self-generation and partly to wave-wave 
turbulent cascading from a large scale, probably corre- 
sponding to the size of SN bubbles. 

Several authors have discussed the possibility that CRs 



could be self-confined by the waves generated through the 
streaming instability that they excite in the direction of 
thei r sp atial gradient (see 0] for reviews) . In particu- 
lar [l3| and [Il| discussed the effect of self-generation in 
the presence of ion-neutral damping and non-linear Lan- 
dau damping (NLLD). The general conclusion in both 
cases is rather interesting: in the Galactic disk and its 
vicinities, waves are damped so fast that CR transport is 
ballistic, while particles are dragged at the Alfven speed 
(confinement) only in the halo, much above or below the 
Galactic disc, where the neutral gas density is sufficiently 
low to avoid effective ion-neutral damping. However, the 
results of these early works were not based on a self- 
consistent solution of the transport equation. More re- 
cently, Ref. [l3] considered the problem of propagation of 
CRs on their way out of the sources (e.g. SNRs): the au- 
thors find that bubbles of self-generated waves develop in 
few hundred pc around SNRs. and that such waves dom- 
inate the diffusive CR propagation. The effective diffu- 
sion coefficient is found to be D{E) oc E'^-^. This slope is 
roughly consistent with the slope of the Boron/Carbon 
ratio, but leads to severe problems with CR anisotropy 
if extrapolated to > TeV energies [ij, [ij]. Moreover 
such approaches cannot explain the complex CR spec- 
trum discussed above since the dependence of D{E) on 
energy does not change. 

Our current knowledge of the Galaxy suggests that, 
although the average density of neutral gas is relatively 
high (~ 1 cm"''), such gas is actually confined to regions 
with small volume filling factor [l^ . Most of the Galaxy 
is filled with a tenuous ionized gas, where ion-neutral 
damping is slow and the wave dynamics is most likely 
determined by NLLD. The implications of the resonant 
absorption of MHD waves by CRs in the Galaxy have 
recently been investigated in [l^ using GALPROP. 

Moreover we have evidence that turbulence exists in 
the ISM with a roughly Kolmogorov-type spectrum up 
to spatial scales of order ^ 50 pc or even larger [l3| . This 
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turbulence is probably injected by SN explosions on such 
scales and then cascades towards smaller scales. It only 
becomes effective for CR scattering when the wavelength 
becomes as small as the particle Larmor radius. Diffusion 
models routinely used in calculations of CR propagation 
are inspired (implicitly or explicitly) by the assumption 
that something like this happens. 

In this Letter we present our calculations of the com- 
bined effect on CR scattering of turbulence cascading 
from some large scale {Lq = 50 pc) through NLLD and 
self-generated waves induced by CR streaming in the 
Galaxy. We find that a change in the scattering prop- 
erties of the ISM must occur at ~ 200 - 300 GV re- 
flecting in a change of shape of the CR spectrum at the 
same rigidity. While the transition energy can be esti- 
mated analytically, we solve the full system of equations 
describing CR transport and wave evolution so as to ob- 
tain a self-consistent spectrum of CRs. In this way, we 
also find that at energies below ^10 GV the advection 
of CRs with waves moving with the Alfven velocity leads 
to a spectral hardening. Both spectral features are ob- 
served [H, HQ : and this work was actually stimulated by 
these observations. 

The calculation — We solve the CR diffusion equation 
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coupled with the equation for the waves: 
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Here f{p, z) is normalized so that the number of particles 
in the range dp around momentum p at the location z is 
Anp^ f{p, z)dp. The diffusion coefficient is related to the 
wave spectrum through the well known expression (lo| : 

where the power in the form of waves, W{k), satisfies: 



Here ^CR is the fraction of the total kinetic energy 
of a SN, ii^sN, assumed to be channelled into CRs, 
and the SN rate is TZsn- The quantity X{ct) — 
47r Jq°° dx x'^~" [Vs;^ + 1 — l] comes from the normaliza- 
tion of the kinetic energy of the SN that goes into CRs. 
Notice that the particle spectrum is assumed to be a 
power law in momentum, as expected for diffusive shock 
acceleration in the test particle regime. 

Eq. (12) describes the stationary wave spectrum W{k) 
under the effect of wave- wave coupling and amplification 
of waves due to streaming instability at a rate Tcrik). 
The cascade is due to NLLD and is described as a diffu- 
sion process in fc-space with a diffusion coefficient [20| : 



Dkk = C^VAk'/^W^kf'^ 



(6) 



for a Kolmogorov phenomenology (Ck ~ 5.2 x 10~^ [13 )■ 
One can easily check that, in the absence of a CR- induced 
contribution, this diffusive process in fc-space leads to 
the standard Kolmogorov spectrum W[k) oc k~^/^ (for 
k S> fco), if the injection of power occurs at a single k^ = 
I/Lq- The effect of CRs is to amplify the waves through 
streaming instability, with the growth rate [lo| : 
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where the spatial gradient of CRs can be found by solving 
the transport equation, Eq. ([1]). 

Eq. ([!]) is solved in the simplifying assumptions that 
D depends weakly on the the z-coordinate, and that the 
Alfven speed is also independent of z, except for the fact 
that Alfven waves move upward (downward) above (be- 
low) the disk. This implies that dvA/dz = 2vaS{z). With 
these assumptions the solution of Eq. ([T|) can readily be 
found to be in the form: 



f{z,p) = fo{p) 
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and fo{p) has to satisfy the following equation, obtained 
by integrating Eq. ([1]) in the range z = (0~ — 0+): 
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with SB'^/Att the power in turbulent fields and Bq the 
regular magnetic field strength. In Eq. ^ the momen- 
tum and wavenumber are related through the resonance 
condition k = ^/^l{p) = 9^0 / {pc) , with the Lar- 
mor radius of particles with momentum p moving in the 
magnetic field Bq. The underlying assumption is that 
SB <C -Bo- In Eq. ([T]) we use a simple injection model in 
which all CRs are produced by SNRs in an infinitely thin 
disc of radius Rd'- 
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The space derivative can be easily derived from Eq. ([8|) 
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A(p) = l-cxp[C(p)]. (10) 



Solving for the CR spectrum in the disc of the Galaxy 
is readily found to be: 
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In the high energy limit, where diffusion prevails upon 
convection at speed va, Eq. (fTTI) reduces to the well 
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known solution of the diffusion equation in one dimen- 
sion, f^'^{p)^qo{p)H/{2D{p)). 

Eqs. (IlJ and ([2]) form a set of two non-linear differential 
equations. Wc solve them itcratively so that the final re- 
sults are the spectrum of CRs in the Galactic disc, fo{p), 
and the power spectrum of waves, W{k), resulting from 
self-generation and cascading. Before illustrating the ex- 
act results, it is useful to estimate the energy of CR pro- 
tons where one expects a transition from self-generated 
waves to waves deriving from Kolmogorov cascade. In 
order to do this, we use the fact that CR driven waves 
saturate when the NLLD rate (Tnl ~ Dkk/k^, with Dkk 
given in Eq. (jH)) equals the growth rate Fcr (Eq. ([7])). Us- 
ing the expression fo{p) = Ap{p/mc)^^p to describe the 
spectrum, with Ap and 7p taken from PAMELA data 
above 250 GeV, we find: 
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where k is related to the particle energy through the res- 
onance condition. The transition length-scale, and hence 
energy, is then found by simply equating W^cR to the ex- 
ternally generated Kolmogorov spectrum. We take the 
latter in the form: 



W^ext(fc) = (2r/B/3fco) (fc/fco)-^/3 



(13) 



Rather than guessing the value of i]b, we use again our 
knowledge of the CR spectrum based on Pamela data, 
and express rys in terms of the CR acceleration efficiency 
^CR- This can be done by just recalling that in the high 
energy regime fa{p) is well approximated by /q '*^(p), with 
D{p) given by Eq. ([3]). At scales at which the Kolmogorov 
turbulence is the dominant scattering source, one finds: 
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Using this condition in Eq. 
Eq. we obtain: 



T3| and equating it to 



228 GeV 



Co. 1 £-5 172-30 



1/3 \ 2(7p-4) 



BoT'\ (15) 



i?d/10kpc, Hs = i?/3kpc, Co.i = 
-BsN/lO^ierg, 7^30 = 72sN/30yr-i, 



where Rd,io = 
Ccr/0.1, E51 = 
Bo,fi = Bo/fiG. 

The estimate obtained for the reference values of the 
parameters is tantalizingly close to the energy where 
PAMELA data show a change of slope of the spectrum of 
protons from i;-2-85±o.oi5 < 230 GeV) to i;-2.67±o.o3 

{E > 230 GeV). Clearly the fact that al low fc's (large mo- 
menta) the power spectrum is W{k) ^ k~^^^ implies that 
the CR injection spectrum must be qo{p) ^ p~^-^ . We 
notice that the transition energy turns out to be indepen- 
dent of the characteristic scale of turbulence, l/k^: this 



fact is especially important given the large uncertainty 
on this parameter. What does depend on feg, linearly, is 
the energy density in turbulent magnetic field, which for 
the above values of the parameters and Lq = 50 pc turns 
out to correspond to « 8%. 

Interestingly enough, taking into account that in the 
advection dominated regime the solution is /o^^(p) ~ 
90 {p) I VA and equating this to the solution in the diffusion 
regime determined by self-generated waves, it is easy to 
see that for reasonable values of the Alfven speed the low 
energy transition to a spectrum foip) oc p~°' occurs at 
E<10 GeV. 

Results and Discussion — The iterative procedure de- 
scribed above leads to the spectrum of CRs plotted in 
Fig. [1] (sohd fine). The symbols are PAMELA data. 
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FIG. 1: Spectrum of CR protons from om calculations (solid 
line) compared with the results from observations of 7-ray 
emission from clouds in the Gould belt [3 (shaded region ) 
and with PAMELA measurements 

Since solar modulation affects the low energy part, only 
the data above 80 GeV were used in [l[ to highlight the 
break at ^ 230 GV. The dashed line is the solution of 
the same equations but setting the self-generation term 
to zero. The shaded area shows the CR spectrum in- 
ferred in Ref. [3. The break at ~ 200 GeV reflects 
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FIG. 2: Diffusion coefficient induced by streaming instability 
of CRs and cascading from a large spatial scale of 50 pc. 

the transition from a regime where the scattering cen- 
ters are self-generated to a regime where particles diffuse 
in external turbulence that cascades from larger spatial 
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scales. This is well visible in Fig.[2l where we plot the dif- 
fusion coefficient with and without self-generation (solid 
and dashed line respectively) . The energy dependence of 
the diffusion coefficient in the energy range 10 < i? < 200 
GeV is not a perfect power law, but if we approximate 
it as such, we find D{E) ~ E'^-'^. At larger energies, 
E > 200 GcV, the trend is D{E) - E^/^, as expected 
for a Kolmogorov cascade. At low energy the change of 
slope is due to the transition to a non-relativistic regime, 
but there the propagation is advection-dominated. 
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FIG. 3: Grammage obtained in our calculations (solid line) 
compared with the leaky box fit proposed in Ref. [10] (dashed 
and dotted lines are the fit for S = 0.5 and 1/3 respectively). 

The grammage traversed by CRs of velocity v 
under the combined effect of diffusion and advec- 
tion at the Alfven speed is easily found to be 
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where fi — 2nmh 



2.4mg/cm is the surface density of the disc (correspond- 
ing to a Galactic disc density Ud — I cm"-^, for a typical 
chemical composition of the interstellar matter). The 
grammage obtained in our calculations is shown in Fig. [3] 
(solid line) and compared with the leaky box fits pro- 
posed in [3 (dashed (dotted) line for S = 0.5 {S = 1/3)). 
The normalization in the 1 — 10 GeV range comes out 
naturally in our calculations, together with the change 
of slope in the diffusion properties at « 200 GeV. 

The flattening in the energy dependence of the diffu- 
sion coefficient at high energies automatically avoids se- 



vere problems with CR anisotropy. The mean anisotropy 
amplitude is (5 ^ 10~^ at 1 TeV and increases with energy 
as E^^^. However, as discussed in [ij], this mean value 
is not very meaningful in that the amplitude is domi- 
nated by the nearest and most recent sources, and may 
dramatically differ from the mean value. Nevertheless, it 
is encouraging that the physical processes described here 
naturally lead to steep energy dependence of the diffu- 
sion coefficient in the low energy regime that does not 
necessarily lead to violate observed data on anisotropy. 

In summary, we showed that in a simplified but physi- 
cally consistent model for the CR propagation in the ISM 
the departure from a power-law spectrum in CRs mea- 
sured at Earth is a consequence of basic processes. Both 
a convective velocity of the order of va in Eq. ([T]) and the 
contribution of CRs to the wave spectrum in Eq. ([5]) are 
unavoidable, albeit often neglected in phenomenological 
studies. It is actually remarkable that without ad hoc free 
parameters, the basic trend shown by the measured and 
inferred CR spectra can be reproduced in a relatively sim- 
ple model. The several approximations that have been 
made (e.g. treatment of non-linear damping, absence of 
re-acceleration, neglect of non-linear effects in diffusive 
shock acceleration) are not expected to change the qual- 
itative picture that emerges. The scenario detailed here 
also allows to accommodate the inferred behaviour of the 
grammage as well as — at least qualitatively — the strin- 
gent constraints coming from anisotropy at high energy. 
We find encouraging that the recent bonanza of CR data 
allows one to gain some insights on (astro)physical pro- 
cesses involving CRs, beyond the mere task of providing 
more accurate fits to injection spectra or propagation pa- 
rameters; this is a trend which hopefully will be further 
boosted in the near future with the expected results of 
AMS-02. 
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